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PREFACE 


The  High  Altitude  Effects  Simulation  (HAES)  Program  sponsored 
by  the  Defense  Nuclear  Agency  since  the  early  1970  time  period,  com- 
promises several  groupings  of  separate,  but  interrelated  technical  ac- 
tivities, e.g.,  ICECAP  (Infrared  Chemistry  Experiments -Coordinated 
Auroral  Program).  Each  of  the  latter  have  the  common  objectives  of 
providing  information  ascertained  as  essential  for  the  development  and 
validation  of  predictive  computer  codes  designed  for  use  with  high  pri- 
ority DoD  radar,  communications,  and  optical  defensive  systems. 

Since  the  inception  of  the  HAES  Program,  significant  achieve- 
ments and  results  have  been  described  in  reports  published  by  DNA, 
participating  service  laboratories,  and  supportive  organizations.  In 
order  to  provide  greater  visibility  for  such  information  and  enhance 
its  timely  applications,  significant  reports  published  since  early  calen- 
dar 1974  shall  be  identified  with  an  assigned  HAES  serial  number  and 
the  appropriate  activity  acronym  (e.g.  , ICECAP)  as  part  of  the  report 
title.  A complete  and  current  bibliography  of  all  HAES  reports  issued 
prior  to  and  subsequent  to  HAES  Report  No.  1,  dated  5 February  1974 
entitled,  "Rocket  Launch  of  an  SWIR  Spectrometer  into  an  Aurora 
(ICECAP  72)",  AFCRL  Environmental  Research  Paper  No.  466,  is 
maintained  and  available  on  request  from  DASIAC,  DoD  Nuclear  Inform- 
ation and  Analysis  Center,  816  State  Street,  Santa  Barbara,  California, 
93102,  Telephone  [805]  965-0551. 


This  report,  which  is  the  first  scientific  report  under  AFCRL 
Contract  No.  F19628-74-C-0177  is  the  fortieth  report  in  the  HAES  series, 
and  covers  technical  activities  period  16  April  1974  through  15  April  1975. 
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The  calculations  in  this  report  provided  a portion  of  the  theoretical 
data  base  on  which  the  EXCEDE  series  of  rocketborne  electron  experiments 
were  designed.  The  support  of  Dr.  Harold  C.  Fltz,  Jr  and  LCDR  Chris 
Thomas  of  DNA  and  Mr.  Herb  Mitchell  of  R&D  Associates  and  the  technical 
direction  of  Dr.  A.T.  Stair,  Jr. , Mr.  James  Ulwick  and  Mr.  Robert  R.  O'Neil 
of  AFGL  is  gratefully  acknowledged. 
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1.0  INTRODUCTION 


Contract  F19628-74-C-0177  is  a contract  to  analyze  High  Altitude 
Effects  Simulation  (HAES)  Experiments  performed  as  part  of  the  ICECAP 
and  EXCEDE  Programs  and  the  application  of  these  measurements  to 
physics  and  chemistry  code  modelling  of  atmospheric  emissions,  such 
as  OPTIR.  The  work  performed  under  the  initial  phase  of  this  program 
includes  modelling  of  energy  deposition,  particularly  for  EXCEDE  type 
experiments  where  radial  distributions  are  probably  more  important  than 
in  auroral  events  (Section  2.  0);  calculations  of  the  zenith  spectral  radi- 
ance at  4.  3 pm  (Section  3.  0);  and  calculations  of  vibrational  population 
of  the  first  excited  level  of  the  CC>2  vibrational  mode  (Section  4.  0). 
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2.  0 


ENERGY  DEPOSITION 


,[1] 


The  energy  deposition  calculations  are  based  on  the  Berger, 


et  alL'J  Monte  Carlo  results  which  are  briefly  summarized  in  the  func- 
tional relations 


F(h,  r)  = ^(h)  R(h.r) 

R(h,  r)  = F(h,  r)/^ F(h,  r’)  2r'  dr' 


A(z  ) = (E/r  ) f(z  /r  ) 
m p m p 


m 


) dh' 


where  F(h,r)  is  the  energy  per  unit  volume  deposited  at  a distance  h-hQ 
down  along  the  field  line  and  r radially  outward  from  a magnetic  field 
line  passing  through  the  injection  point  of  a single  electron.  Needed  to 
calculate  this  are  h),  the  atmospheric  density  at  the  point  of  observa- 
tion; R(h,  r)  the  radial  distribution  function;  z , the  mass  thickness  be- 

m 

tween  the  points  of  injection  and  observation;  and  r , the  practical  range 

P 

of  the  incident  electron  (Figures  1,  2,  3). 

The  Monte  Carlo  calculations  have  been  fitted  by  Berger,  et  al, 
to  the  following  functional  approximations: 

a-2 

J *lL  f_r\  e-b(r/r„) 

rH  ' 


R(h,r) 


2*r 


f(  a) 


H 


H 


where 


r - the  initial  Larmor  radius  of  the  incident  electron 
H 

a,b  - dimensionless  parameters  related  to  moments  of  R(h,  r) 

2 2 2 2 

where<^3>is  the  mean  radius  and  the  variance  is  (f;  cr  -<p  > ~ <p> 

<Q>2  ■ 

a * 0-2  ’ b - <p> 

Figure  4 is  a plot  of  1/Va  as  a function  of  z /r  . Figure  5 is  a plot  of  a/b 


m 

as  a function  of  energy  for  different  values  of  z /r  . 

° m p 

of  a/b  as  a function  of  z /r  for  3 keV  electrons. 

m p 


Figure  6 is  a plot 
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Tables  in  the  Berger,  et  al  paper  provide:  radial  distribution 

a and  b parameters  for  a range  of  energies  and  distances  from  the  point 

of  injection:  atmospheric  density  as  a function  of  altitude;  and  values  ot 

the  "universal"  f-function  of  the  axial  factor,  A(z  ),  for  various  dis- 
til < 

tances  from  the  source.  These  tabular  entries  will  be  least- squares 
fitted  to  appropriate  functions  and  incorporated  into  a subroutine  that  re- 
turns the  energy  deposition  rate  and  ionization  production  rate,  Q = n^  = 

F(h,  r)  • I /ew  (I  - beam  current,  e - electron  charge,  w - eV/ion  pair), 
o o 

for  a given  energy  and  geometry. 

Recombination 

For  sufficiently  high  ionization  densities,  recombination  may 
modify  results  significantly  and  thus  must  be  taken  into  account. 

The  basic  equation  relating  the  ionization  density  to  source  strength 
is: 

h (t)  = Q(t)  - dn2(t) 
e e 

a nonlinear,  first-order  differential  equation  of  the  Riccatti  type.  The 
analytical  solution  for  arbitrary  Q(t)  is  not  known,  but  through  the  non- 
linear transformation 

ft(t)  = exp  (a  \ dt'  n (t)) 
e 

The  linear,  but  second-order,  homogeneous  differential  equation 

A 

n (t)  - CtQ(t)  ft  (t)  = 0 
e e 

is  obtained,  which  for  certain  source  term  functional  dependences,  can 
be  solved  in  terms  of  known  functions. 

Line  Integrals 

The  electron  density,  n^fr,  t),  being  defined  by  a nonlinear  dif- 
ferential equation  involving  a time -varying  source  function,  will  virtually 
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never  be  obtained  in  an  analytical  form  so  that  a numerical  evaluation  of 

the  line  integral  fn  (r,  t)dl,  is  unavoidable.  For  each  sub-interval  of  the 
J e 

integral,  the  value  of  n (r,t)  will  have  to  be  obtained  by  an  integration  of 

e 2 

the  equation  n^t)  = Q(t)  -anjt),  where  the  time  variation,  Q(t),  of  the 
source  is  due  to  source  trajectory  motion  and  possibly  to  source  inten- 
sity changes. 

Approximate  Calculations 


For  order  of  magnitude  and  qualitative  investigations,  certain 
simplified  cases  have  been  studied.  If  recombination  during  the  rela- 
tively brief  duration  of  the  source  pulse,  Q(t),  can  be  neglected,  then 


n (t)  = Q(t)  - an  (t)  ~ Q(t) 
e e — 


and  the  electron  density  during  the  source  pulse  is  given  by 


ne(t)  -/ dt'  Q(t,) 


After  the  source  pulse  is  over,  n = -an  , which  has  the  solution 

e e 


e ' 1 +an  t 

o 

-CO 

where  n = / dt 1 Q(t')  is  the  total  electron  density  produced  by  the  sourc 

° TOO 

pulse  (in  the  absence  of  recombination).  Since,  for  the  general  point  in 
space, 


Q(t)  = Q(h(t),  r(t))  = I F(h,  r) / ew 

o 

where  h and  r are  determined  by  the  (moving)  source  and  observation 
point  geometry,  then 

n (?) 

, -*■  , , o 

n (r , t)  = 7—7 

e 1 + an  (r)t 

o 
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where 


n (r)  = -2-  J dt'  F(h(t),  r(t)) 
o ew  */qd 

To  perform  an  analytical  line  integration  of  this  n (? , t),  the  "initial" 

e 

distribution,  n^(r),  must  be  of  relatively  simple  algebraic  form.  Un- 
fortunately, the  n (?)  from  the  Berger  theory  is  generally  a non-element- 
ary  integral,  effectively  making  a closed  solution  impossible.  However, 
in  the  special  case  of  line  integrals  along  directions  normal  to  the  plane 
of  the  magnetic  field  and  source  velocity,  and  at  altitudes  for  which  the 
radial  distribution,  R(h,  r),  has  an  integer  a-parameter,  nQ(r)  can  be 
evaluated  in  terms  of  the  zero  and  first-order  modified  Hankel  functions. 
These  special  n (r)  cases  can  be  evaluated  from  tables  and  least- squares 
fit  to  functions  that  permit  analytical  evaluation  of  the  line  integral  to 
yield  the  approximate  time  dependence  of  Jn  dl  at  the  special  altitudes. 

The  model  distributions  which  were  fit  to  the  special  (evaluable) 
integrals  were  rectangular  (uniform),  triangular,  exponential  and  Lor- 
entzian.  The  resulting  time  profiles  differed  by  40%  at  one  characteris- 
tic time,  T ; by  3:1  at  10  f and  by  10:1  at  100  X • 
o o'  o 

More  Exact  Calculations 

With  some  idea  of  the  source  function  pulses,  Q(t),  to  be  en- 
countered, the  electron  density  n^(r,  t)  can  be  estimated,  now  taking 
recombination  into  account,  either  by  numerically  integrating  the  basic 
equation 

n (t)  = Q(t)  - an2  (t) 
e e 

defining  Q(t)  directly  in  terms  of  the  Berger  function,  F(h,  r),  or  by 
segmental  representation  of  Q(t)  with  polynominal  fits  to  selected  values 
of  the  Berger  function. 

There  exist  a number  of  straightforward,  first-order  differential 
equation  numerical  integration  techniques,  but  good  numerical  accuracy 
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t 


1 


requires  evaluation  of  a relatively  complicated  Q(t)  function  many  times 
over.  Since  the  many  such  time -integrations  required  in  the  various  de- 
sired line  integrals  might  prove  too  costly  in  computer  time,  a linear  seg- 
ment representation  approach  was  also  looked  into.  If  the  Q(t)  function 
can  be  adequately  represented  by  five  or  ten  linear  segments,  then  the 
solution  throughout  any  segment  can  be  expressed  in  terms  of  the  solu- 
tions of  the  transformed,  second-order  differential  equation: 

n (t)  + C(  Q(t)  n (t)  = 0 
e e 

where 

Q(t)  = Q (1  + t/t) 
o 

the  solutions  of  which  are  directly  related  to  the  Airy  integrals  Ai(x) 
and  Bi(x)  and  their  derivatives.  Subroutines  to  evaluate  these  functions 
for  -oo<x<oo  have  been  written  and  await  testing. 

The  energy  deposition  calculations  discussed  above  are  applicable 
to  both  the  ICECAP  and  EXCEDE  Programs  and  will  be  used  to  evaluate 
electron  spectra  such  as  that  presented  in  Appendix  1.  The  accuracy  of 
the  radial  description  is  more  relevant  to  the  EXCEDE  geometry  than  to 
a typical  auroral  geometry. 


| 
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3.0  SPECTRAL  RADIANCE  CALCULATIONS  AT  4.  3 UM,  9.  6 JJM, 

! 15  UM. 

The  night  time  zenith  spectral  radiance  has  been  calculated  for  the 
CO2  bands  at  4.  3 pm  and  15  pm  and  for  the  ozone  band  at  9.  6 pm  at  alti- 
tudes for  60  to  120  pm.  The  spectral  radiance  was  calculated  assuming 
the  U.  S.  Standard  Atmosphere  (180°  K Mesopheric  temperature)  and  an 

atmosphere  with  a 230°  K Mesopheric  temperature.  The  results  are  pre- 

2 

sented  in  Appendix  2.  Degges  1974  Infrared  Radiance  Model  was  used  for 
these  calculations. 

, 3 

Degges'  work  is  an  extension  of  the  study  of  Corbin,  et  al.  (1969) 
and  Degges  (1972)?  The  former  investigated  the  natural  infrared  back- 
ground of  the  earth  in  the  5 to  25  micrometer  spectral  region.  For  con- 
venience, their  study  divided  the  atmosphere  into  two  regions  with  a divi- 
sion at  70  km.  Below  70  km  the  atmosphere  was  assumed  to  be  in  thermal 
equilibrium.  Above  70  km  explicit  calculations  were  made  of  processes 
which  excite  and  de-excite  molecular  vibrational  and  rotational  levels  which 
are  the  source  of  infrared  radiation.  Their  study  concentrated  on  radiation 
from  water  vapor,  carbon  dioxide,  ozone,  nitric  oxide  and  nitrous  oxide, 
which  are  the  principal  radiating  species  in  the  spectral  region  considered. 

In  addition,  nitric  acid  was  included  in  the  lower  atmosphere  work  and 
estimates  were  made  of  radiation  to  be  expected  from  particulate  matter 
suspended  in  the  atmosphere. 

Corbin,  et  al.  (1969)  presented  models  for  the  lower  atmosphere 
for  a wide  range  of  seasonal  and  latitudinal  conditions.  This  was  not 
possible  for  the  abundances  of  most  minor  neutral  species.  More  data  has 
since  become  available,  particularly  for  nitric  oxide  and  the  hydroxyl 

radical,  but  at  present  it  appears  that  the  best  means  of  estimating  abun-  ! 

dances  of  important  infrared  emitting  species  is  chemical  rate  equation 
integrations  including  molecular  diffusion  and  eddy  mixing.  Even  calcula- 
tions involving  transport  properties  cannot  always  be  accepted  because  the 
values  of  eddy  mixing  coefficients  are  to  a large  extent  only  informed 
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guesses  and  published  calculations  often  use  outdated  rate  coefficients. 
Degges  (1972)  reported  a computer  program  with  which  to  determine 
diurnal  variations  in  abundances  of  minor  species,  to  investigate  the 
effects  of  changes  in  assumed  eddy  mixing  coefficients,  and  to  estimate 
the  effects  of  new  determinations  of  chemical  rate  coefficients. 


A second  area  of  study  involved  in  improving  the  radiance  model 
includes  the  physical  processes  that  control  the  population  of  infrared 
emitting  states  of  atmospheric  molecules.  Except  for  the  pure  rotational 
radiation  from  molecules  such  as  water,  the  degree  of  excitation  of 
vibrational  levels  determines  the  radiation  from  infrared  emitting  mole- 
cules. The  most  important  mechanisms  are  collisional  excitation  and 
de-excitation  and  absorption  and  re-emission  of  electromagnetic  radiation. 

In  the  troposphere  and  lower  stratosphere,  collisional  processes  are 
rapid  enough  to  control  the  population  of  vibrational  levels.  Above  30  to 
50  km,  however,  collisional  excitation  becomes  less  efficient  and  radia- 
tive processes  become  important.  The  combined  effects  of  collisional 
and  radiative  processes  must  therefore  be  considered.  Below  an  altitude 
of  about  90  km  molecular  nitrogen  and  oxygen  are  the  most  important 
collision  partners.  Above  that  altitude,  atomic  oxygen  becomes  impor- 
tant, both  in  exciting  nitric  oxide  and  in  determining  molecular  oxygen  and 
nitrogen  vibrational  temperatures. 

In  determining  the  effects  of  radiation  on  the  populations  of  vibra- 
tional levels  of  infrared  active  molecules,  it  is  necessary  to  separate  the 
radiation  of  a single  change  in  vibrational  quantum  numbers  from  the  rest 
of  the  radiation  field.  The  previau  sly  reported  studies  did  this  by  assum- 
ing a Doppler  line  shape  for  the  individual  rotational  lines  of  a band  and 
were  able  to  obtain  adequate  numerical  approximations  for  radiative  trans- 
fer functions  appropriate  to  single  bands  of  linear  molecules,  and  less 
accurately,  for  water  vapor  and  ozone  bands  lying  in  the  spectral  region 
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4.0  MODEL  CALCULATIONS  OF  THE  VIBRATIONAL  POPULATION  OF 


THE  FIRST  EXCITED  LEVEL  VIBRATIONAL  MODEL. 

Ln  order  to  gain  insight  as  to  possible  effects  of  the  time  and  space 
variations  in  auroral  excitation  of  the  carbon  dioxide  (CO^)  infrared  band 
systems,  we  have  continued  model  calculations  of  the  vibrational  popula- 
tion of  the  first  excited  level  of  the  CO2  V3  vibrational  mode.  The  model 
atmosphere  used  contains  only  CO2,  molecular  nitrogen,  and  atomic  and 
molecular  oxygen.  The  C mixing  ratio  is  assumed  constant  at  320 
parts  per  million,  by  volume.  Calculations  are  restricted  to  levels  two 
kilometers  apart  lying  between  50  and  150  kilometers  altitude.  The  tem- 
perature falls  from  270  K at  50  km  to  a mesopause  temperature  of  191 
between  80  and  90  km.  Above  90  km,  the  temperature  rises  to  66 1 K at 
1 50  km. 

The  radiative  transfer  model  of  Degges  (1974)  is  used  with  the 
collision  processes  of  Table  1,  the  rates  shown  being  those  given  by 

5 

Garvin  and  Hampson  (1974).  The  QCHEM  integration  algorithm  described 
by  Manley,  et  a 1 (1973)  is  used  to  integrate  rate  equations  for  the  vibra- 
tional populations  of  the  CO2  mode  and  the  first  vibrational  levels  of 
N£  and  C>2.  The  results  presented  in  the  following  figures  should  be 
regarded  as  only  semi-quantitative.  The  system  of  rate  equations  is  stiff 
with  coupling  between  stiff  components  through  the  radiation  field.  Manley, 
et  al  (1973)  point  out  that  in  cases  of  this  nature,  the  solution  given  by  the 
QCHEM  algorithm  relaxes  in  the  correct  direction,  but  the  numerical 
results  are  not  properly  correlated  with  the  prescribed  time  steps. 

Results  are  presented  here  for  two  (2)  sample  cases.  In  the  first, 
vibrational  temperatures  at  all  altitudes  were  set  to  250  K initially. 
Assuming  night  time  conditions,  the  system  was  integrated  for  10^  seconds. 
Then,  an  auroral  "drizzle"  was  introduced  with  a peak  ionization  rate  of 
2.  5 x 10^  ion  pairs  cm'  sec  at  110  km  and  the  integration  continued 
for  10^  seconds.  Then,  a strong  arc  was  turned  on,  with  a peak  ionization 
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rate  of  10^  ion  pairs  cm"  5 sec  " * at  100  km  and  the  integration  continued 
for  10^  seconds.  Following  this,  the  auroral  excitation  was  turned  off  and 
relaxation  followed  for  10^  seconds.  The  procedure  was  repeated  for  an 
initial  vibrational  temperature  of  300  K.  In  both  cases,  the  assumption  was 
made  that  for  each  ion  pair  produced,  three  nitrogen  molecules  are  excited 
to  the  first  vibrational  level. 

Figures  1 through  4 display  results  for  the  case  where  initial  vibra- 
tional temperatures  are  250  K.  Figures  5 through  8 display  results  for  the 
case  where  initial  vibrational  temperatures  are  300  K. 

In  Figure  1,  vibrational  populations  increase  with  time  below  60  km 
as  collisions  and  radiative  transport  are  effective  in  bringing  vibrational 
temperatures  closer  to  kinetic  temperatures.  Above  75  km,  populations 
fall  with  time.  At  late  times,  populations  rise  above  125  km  as  collisions 
are  effective  in  increasing  nitrogen  vibrational  temperatures  toward  kine- 
tic temperatures,  with  subsequent  transfer  of  vibrational  excitation  to  C02. 

There  is  little  change  in  Figure  2,  the  "drizzle"  excitation  not  being 
large  compared  to  the  normal  collisional  and  radiative  mechanisms. 

The  auroral  effects  are  clearly  evident  in  Figure  3.  Above  100  km,  a 
balance  is  reached  between  enhanced  auroral  excitation  and  radiative  decay 
by  10^  seconds.  This  aurora  produces  no  ionization  at  90  km  and  below. 

At  late  times,  vibrational  populations  continue  to  increase  at  altitudes  as 
low  as  85  km  as  energy  is  transported  down  by  radiation.  This  downward 
transport  of  radiation  continues  after  the  aurora  is  turned  off,  as  shown  in 
Figure  4. 

The  behavior  for  the  higher  initial  vibrational  temperature  of  300  K 
is  similar,  as  shown  in  Figures  5 through  8,  except  that  vibrational  tem- 
peratures fall  at  all  times  for  altitudes  below  65  km.  The  auroral  effects 
are  not  as  visibly  pronounced,  because  of  the  higher  initial  vibrational 
temperature. 

Figures  9 through  12  present  vibrational  temperatures  corresponding 
to  the  vibrational  populations  of  Figures  1 through  4. 
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COz  vibrational  populations  at  10  ,10  .....  10  , 10 

seconds  following  relaxation  from  a uniform  vibrational 
temperature  of  250  K. 

CO^  vibrational  populations  at  10  ^ , 10  ^ , . . . . 10"\  10^ 
seconds  with  auroral  "drizzle",  following  conditions  of  Figure  I. 
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CO^  vibrational  populations  at  10  ,10  ,....10,10 

seconds  with  strong  auroral  arc,  following  conditions  of  Figure  2 
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CO^  y ^ vibrational  populations  at  10  , 10  ,....10  , 10 

seconds  during  relaxation  from  conditions  of  Figure  3. 
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CC>2  vibrational  populations  at  10  ,10  ,....10,10 

seconds  following  relaxation  from  a uniform  vibrational 
temperature  of  300  K. 
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CO^  V3  vibrational  populations  at  10  ,10  ....  10  , 10 

seconds  with  auroral  "drizzle",  following  conditions  of  Figure  5. 
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CO^  V3  vibrational  populations  at  10  ,10  ,....10,10 

seconds  with  strong  auroral  arc,  following  conditions  of  Figure  6 
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COz  vibrational  populations  at  10  ,10  ,....10  , 10 

seconds  during  relaxation  from  conditions  of  Figure  7. 


CC>2  vibrational  temperatures  corresponding  to  vibrational 
populations  of  Figure  1. 

CO,  vibrational  temperatures  corresponding  to  Figure  2. 
C02  y vibrational  temperatures  corresponding  to  Figure  3. 
CC>2  V'3  vibrational  temperatures  corresponding  to  Figure  4. 
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Appendix  1 


Auroral  Electron  Energy  Spectra 
Rocket  A18.  219-1 
Launched  25  February  1974 
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The  data  in  the  following  figures  was  presented  at  the  DNA  ICECAP 
Data  Meeting  in  January  1975,  Rocket  A 18.219-1  was  launched  25  Febru- 
ary 1975  and  overflew  the  same  bright  auroral  arc  twice.  Peak  energy 

-2  - 1 

deposition  was  in  excess  of  100  ergs  cm  sec  . The  electron  energy 
deposition  was  measured  with  two  instruments,  an  aluminum  covered 
scintillator  coupled  to  a photomultiplier  and  an  electrostatic  analyzer 
(ESA).  Figure  1 shows  schematically  the  operation  of  the  ESA.  Spherical 
analyzer  plates  are  used  to  select  a narrow  energy  band.  The  electrons 
are  then  accelerated  with  a voltage  applied  to  a thin  aluminum  deposition 
on  the  face  of  a calcium  fluoride  scintillator.  This  extends  the  low  energy 
response  of  the  instrument.  By  modulating  the  analyzer  plate  voltage 
and  using  a narrow  band  amplifier  the  sensitivity  of  the  instrument  is 
significantly  improved.  Figure  2 shows  the  scintillator  measurement 
as  a function  of  time.  This  detector  measures  the  integrated  electron 
flux  above  4 or  5 keV.  The  two  peaks  correspond  to  the  time  when  the 
rocket  is  on  the  magnetic  field  line  connecting  to  the  visible  aurora.  The 
auroral  motion  is  presently  being  analyzed  as  part  of  the  data  analysis. 

The  visible  aurora  initially  is  to  the  north,  moves  to  the  south  under  the 
rocket,  and  then  moves  to  the  north  again  passing  under  the  rocket. 

The  following  figures  show  the  auroral  energy  spectrum  measured 
with  the  ESA.  The  insert  shows  the  scintillator  measurement  and  the  in- 
verted triangle  indicates  when  that  electron  energy  spectrum  was  taken. 

The  auroral  peak  energy  extended  beyond  30  keV  and  showed  significant 
spectral  differences  north  and  south  of  the  peak  region.  The  spectral 
form  north  of  the  peak  region  was  much  softer  than  that  south  of  the  peak 
region.  These  spectral  differences  persisted  throughout  the  flight.  The 
most  striking  feature  is  that  the  electron  spectral  description  of  the  auroral 
form  varied  very  little  while  the  form  moved  rapidly  in  the  sky.  The  de- 
tailed analysis  of  this  aurora  will  be  presented  in  a later  report. 
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Appendix  2 

Spectral  Radiance  Calculations  at 


4.  3 pm,  9.  6 pm,  and  15  pm 
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Figure  A.  2-1  Spectral  Radiance  in  Upward  Direction 
above  indicated  altitude,  US  Standard 
Atmosphere,  1962.  CO^  4.  3 um  band. 
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Figure  A.  2-2  Spectral  Radiance  in  Upward  Direction 


above  indicated  altitude.  230  K Mesopause, 


Figure  A. 2- 3 


Spectral  Radiance  in  Upward  Direction 
above  indicated  altitude,  US  Standard 
Atmosphere,  19&2.  CO^  3 ^ m ^and. 
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Figure  A.  2-4  Spectral  Radiance  in  Upward  Direction 

above  indicated  altitude.  230  K Mesopause, 
4.  3 u m CO^  band. 
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Spectral  Radiance  in  Upward  Direction 
above  indicated  altitude,  US  Standard 


Atmosphere,  1962.  Ozone  9.6  u m band 
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